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We introduce the Bioj3 Framework, a meta-model for both protein-level and membrane-level interac- 
tions of living cells. This formalism aims to provide a formal setting where to encode, compare and 
merge models at different abstraction levels; in particular, higher-level (e.g. membrane) activities can 
be given a formal biological justification in terms of low-level (i.e., protein) interactions. 

A Bioj3 specification provides a protein signature together a set of protein reactions, in the spirit 
of the K"-calculus. Moreover, the specification describes when a protein configuration triggers one of 
the only two membrane interaction allowed, that is "pinch" and "fuse". 

In this paper we define the syntax and semantics of Bioj3, analyse its properties, give it an in- 
terpretation as biobigraphical reactive systems, and discuss its expressivity by comparing with K- 
calculus and modelling significant examples. 

Notably, Bioj3 has been designed after a bigraphical metamodel for the same purposes. Hence, 
each instance of the calculus corresponds to a bigraphical reactive system, and vice versa (almost). 
Therefore, we can inherith the rich theory of bigraphs, such as the automatic construction of labelled 
transition systems and behavioural congruences. 

1 Introduction 

Cardelli in [5] has convincingly argued that the various biochemical toolkits identified by biologists can 
be described as a hierarchy of abstract machines, each of which can be modelled using methods and 
techniques from concurrency theory. Like other complex situations, it seems unlikely to find a single 
notation covering all aspects of a whole organism; instead, several models, often presented as process 
calculi, have been proposed in literature, each focusing on specific aspects of the biological system, at 
different levels of abstractions. This arises the problem of how to integrate these models. Indeed, these 
machines operate in concert and are highly interdependent: "to understand the functioning of a cell, one 
must understand also how the various machines interact" [5]. To this end, we need a general metamodel, 
that is, a framework, where these models (possibly at different abstraction levels) can be encoded, and 
their interactions can be formally described. Eventually "we will need a single notation in which to 
describe all machines, so that a whole organism can be described" [5]. 

As a step towards this long-term goal, in this paper we present Bio/3, a framework calculus for 
dealing with both protein-level and membrane-level interactions. Many specific models, with protein 
interactions and membrane trafficking, can be encoded in Bio/3 by providing a specification. Notably, 
different specifications can be merged, which corresponds to "put together" different models and systems, 
often at different levels of abstraction. This allows to investigate the interactions between these models, 
which otherwise would be kept separate; for instance, this allows to foresee emerging properties, such 
as behaviours due to interactions of different abstraction levels, and which cannot be observed within a 
single machine model due to its intrinsic abstractions. 

A fundamental design choice is that this framework has to be biologically sound, i.e., it must admit 
only systems and reactions which are biologically meaningful, especially at lower level machines (i.e. 
protein). This is important because in this way, encoding a given model in the framework provides 
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automatically a formal, biologically sound justification for the model (or "implementation") in terms of 
protein reactions and explains how its membrane-level interactions are realised by protein machinery. 
Also, the representation of a specific model is less error-prone. 

To this end, Bio/3 has been designed after a bigraphical metamodel, called biobigraphs and biological 
bigraphical reactive systems (BioRS) and presented in a companion paper [3], for dealing with both 
protein-level and membrane-level interactions. Actually, Bio/3 can be seen as a syntactic formalism for 
representing precisely biobigraphical reactive systems. However, for sake of simplicity and lack of space, 
in this paper we will not discuss the connection between Bio/3 and biobigraphs, which will be presented 
in forthcoming work. For the purposes of the present paper, it suffices to know that biobigraphs, and 
hence Bio/3, are adequate with respect to Danos and Laneve's K-calculus, one of the most accepted 
formal model of protein systems: we can describe all and only protein configurations and interactions of 
the fc-calculus. (Is important to notice, however, that our methodology is general, and can be applied to 
other formal protein models.) Hence, the Bio/3 Framework can be seen somehow as an extension of the 
fc-calculus, adding biological compartments (that is membranes); however, this is a consequence of the 
fact that the underlying model, the biobigraphs, are adequate w.r.t. the fc-calculus. 

When membranes come into play one wants to deal also with membrane transport, hence, we have to 
take into consideration reactions like endocytosis, exocytosis, cellular fusion and vesiculation. Actually, 
as observed by biologists (see [1] and [2, Ch. 15]), membrane interactions present always a "preparation 
phase", where some proteins (receptors and ligands) interact to get in place, followed by the actual 
"membrane reconfiguration" phase, which can be either a pinch or a fuse [6]. Therefore, biobigraphs, 
and Bio/3, provide just two general rules for all membrane interactions, whilst the preparation phase 
depends on the specific proteins involved and hence left to the encoding of the specific model. 

Summarizing, for encoding a given model one has just to instantiate this general schema by specify- 
ing (a) a protein signature, that is a set of abstract proteins; (b) a set of protein rules, describing protein 
interactions ignoring compartments; (c) a set of protein configurations which trigger a mobility reaction. 
Then, all other reactions (mobility, in particular) are automatically provided by the framework, in a sound 
way. Hence, modelers can focus on protein-level aspects, leaving to the framework the burden to deal 
with biological compartments and membrane transport. 

A further motivation for our framework comes from the many general results provided by bigraphs. 
We mention here only the construction of compositional bisimilarities [10], allowing to prove that two 
systems are observational equivalent, that is, they can be exchanged in any organism without that the 
overall behaviour will change. Also for this application it is important to restrict the systems allowed by 
the framework to only those biologically meaningful. 

Synopsis In Section 2 we present the Bio/3 framework: its syntax, a type system for characterizing 
well-formed (i.e., biologically meaningful) terms, and general operational semantics. As an example 
application of this framework, in Section 3 we model the mechanism of membrane traffic. In Section 4 
we provide a formal connection between fc-calculus and Bio/3. Related work is discussed in Section 5, 
and in Section 6 we draw some conclusion and plan future work. 

2 The Bioj3 Framework 

In this section we present the Bio/3 framework. Although one should not think of Bio/3 necessarily as an 
extension of the fc-calculus, its connections with that model are quite strong, and hence it is convenient 
to reuse notations and concepts from K. 
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vesciculation endocitosys fusion exocitosys 



Figure 1: Vesiculation and endocytosis can be generalized by means of an inward and outward pinch, 
and fusion and exocytosis by a horizontal and vertical fuse, respectively. 

2.1 Syntax 

The language of the Bio/3 framework is parametric over a protein signature (3? = S? p U 2? a p-,s '■ & — ► 
N), where is a finite set of (abstract) proteins, partitioned into polar £P p and polar 3P ap . JV is an 
enumerable set of edge names, and s : & — ► N assigns to every protein the number of its domain sites. 

A (protein) interface is a map from N to JV + {h, v} (ranged over p, a, .. .). Given an interface p 
and a protein name A G a site (A,i) is visible if p(i') = v, hidden if p(z) = /z, and fz'ecf if p(z) G ,/f . 
A site is free if it is visible or hidden. In the following, we will write p = 1 + 2 + 3* to mean p(l) = v, 
p(2) = h, and p(3) = x. Denote with \p,x\ the number of occurrences of x in p, that is, \p,x\ = |{z G 
s(p)\p(i)=x}\. 

The syntax of Bio/3 over a given protein signature, consists of systems, which can be nested, and 
membranes. Polar proteins can freely float at the level of systems, while apolar proteins can be embedded 
into membranes. Let A t G & t , for t G {p, ap}, n G , and p : N — » + {/i, v} be a protein interface. 

P,e::=o|A p (p)|^lPJ5|P*e|vn.P|p„?P|f n ^5lPjj (Systems) 
5, T ::= | A ajJ (p) | S*r | pj- ;5 | (Membranes) 

Systems are polar solutions (soluble in aqueous environments), which can be either the empty system 
o, or a polar protein A p (p ) , or a compartment IS \ Pjj, or a group of systems P * Q, or a system prefixed 
by a new name constructor vn.P. Membranes are apolar solutions (soluble in oily environments), which 
can be the empty membrane 0, or an apolar protein A ap (p), or a group of apolar solutions S* T . Notice 
that membranes are not sequences of actions (as, e.g., in Brane calculus [4]) but simply a collection of 
proteins. As usual the "new" operator is a binder: in vn.P, P is the scope of the binder vn. Sharing of 
names on interfaces represent protein domain-domain bonds. For instance (x)(A(l x ) ,S(1 +2 X ))) is a 
protein solution where there is a bond between sites (A, 1) and (B, 2). 

Finally, the special constructors pinch p and fuse f, together with their co-actions p ,f , to represent 
when a compartment reconfiguration takes place. These actions are inspired by the two basic membrane 
reconfiguration processes, namely "pinch" and "fuse". A membrane can be pinched (either outward or 
inward, Figure 1 left), forming a bubble that detaches producing a new compartment; or two membranes 
can fuse, in either directions (Figure 1 right) mixing contents. Intuitively, p„ , P can be seen as a system 
P waiting to complete, or commit, a pinch interaction with another system featuring the co-action p^-; 
similarly for f . Thus, these actions are prefixes which freeze the continuation system/membrane until the 
compartment reconfiguration is completed. The subscripts pair-up corresponding actions and co-actions. 

Let us define the co-action operator mapping each mobility action to its dual, and such that 
((•)- L )- L = id; this operator is extended to sets of actions as S 1 - = {s 1 - \ s G S}. Moreover, for a set 
of names X C ,jV , the restriction of a set of actions over X is S\ x — {s x € S | x G X}. In the following, 
P, Q,R, ■ ■ ■ range over systems, S,T,... over membranes and K,L,H,... over both systems or membranes. 
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The set of free names fn(P) and fn(S) are defined as usual, covering also subscripted names on 
mobility actions (clearly restriction on names binds also action names). In the following fn(K\ , . . .,K n ) 
stands for (J" =1 fn(Ki), for every K{ either a system or a membrane. 

Let Act = { 

Pn j P« ) f/i > fn" I w ^ ^} be the set of all mobility actions. We define the sets of occurring 
actions act(P), act{S) C ACT on systems and membranes, respectively, as follows: 

acf(o) =at*(Ap(p)) = act{lS\P]]) = act (S) Li act (P) act(P*Q) = act(P)L)act(Q) 
act(0) = act{A ap {p)) = actiyn.P) = act(P) act(S*T) = act(S) Li act (T) 
act(p n ° 9 P) = {p n }L)act(P) act(f„$P) = {f n }U act (P) act(p^° 9 S) = {p,f } Dact(S) act(f„) = {f^} 

Bio/3 terms can be rearranged according to a structural equivalence =, defined to be the least equiv- 
alence closed under a -equivalence and satisfying the rules below: 

OP € {*,*} OP £ {*,*} 

P*o = P 5*0 = 5 KopL = LopK Kop (Lop//) = (K op L) opH 

nj-m n j fn{Q) n $ fn(S) 

vn.vm.P = vm.vn.P Vn.o = o vn.(P* Q) = (vn.P) * Q Vn.lS\P^ = lSlvn.P^ 

opG{*,*} K = L P = Q S = T ?G{p,pV} K = L ?G{p,f} n^m 

KopH = LopH IS\P^ = ITIQ^ t n %K = t n %L Vn.(t m %P) = t m %Vn.P 

Well-formedness As for other languages for graph-like structures, the syntax above is too general, 
since many syntactically correct terms do not have a clear biological meaning. In this subsection give 
an informal description of the well-formedness conditions that terms must satisfy; in the next subsection 
we will present a type system enforcing these conditions. 

Unsurprisingly, these conditions are similar to those of related models, such as the fc-calculus, since 
they are motivated by biochemical properties of proteins and membranes. 

First, protein bonds can form only between exactly two active sites on protein surfaces. Hence, terms 
like (x)(A(l x )) or (x)(A(l x ) ,A(1 X ) ,B(l x )) have to be discarded. This requirement is extended to action 
names, because we want to pair-up just a couple of actions at once. Furthermore, a pinch action can be 
connected only to co-pinch, and similarly for fuse and co-fuse 1 . 

Regarding membranes, we have to ensure the impermeability of membrane boundaries: complexa- 
tions cannot take place through membranes, hence protein bonds cannot cross them. Thus, terms like 
A( V) * £0 1 B( have to be ruled out. 

Another (more technical) condition is that continuations after mobility actions cannot have other 
actions. The intuition is that a membrane engaged with a pinch or fuse interaction, has to complete it 
before initiating another one. 

We summarize these conditions in the follwing definition. 

Definition 2.1 (Well-formedness). A system P is well-formed if the following conditions hold: 

Graph-likeness: free names occurs at most twice in P, and every binder ties either or 2 occurrences; 
Impermeability: for any occurrence of a compartment \S\ Q]] in P, fn(Q) C fn(S) hold; 
Action pairing: for any occurrence of an action t n (t G {p,f}), the name n can appear only in the co- 
action t^- 



1 Actually, it suffice to disallow name sharing between actions and protein interfaces, but strengthening the condition will 
become useful in proving the subject reduction theorem. 
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. , ££{0,0} Ae^ Vxefn(p). \p,x\ <2 , , 

( empty) L J v ; 1 — (prot) 

1 Hyj 0;0H £ :0 {x G /«(p) | |p,x| = 1};{x G /ra(p) | \p,x\ = 2} H A(p) : ; 

, . , fe{p, P \f} r i; r 2 H^:0 ac*(£) = r^hP:? x^n t\ {x} = 9 
(actlon) f^T^-iZTfe} r i: r 2 \Whv^:r (v " prot) 

(co-f) , r- ^{ P ,f} r i; r 2 ,*hP:TU{^} {^}nr = e 

V x;0hf^:{f^} r i; r 2 hvx.P:T v ; 

opG{*,*} T[,T;T2\~ K : T Ai,r;A 2 h L : o T u T\T 2 'r S : z T;A 2 hP:a 

ok e (Ti,r 2 ,Ai,A 2 ) ok r (r,T,a) ok e (ri,r 2 ,0,A 2 ) ok r (r,T,a) 

par cell 

VH 7 ri,Ai;r 2 ,A 2 ,rh^opL:TUa ri;r 2 ,A 2 ,ri- IS\PN : TU<7 v 7 



Figure 2: Type system for Bio/3 terms. 



Action prefix: for any occurrence of a prefix , 2 (t G {p, p^, f}) in P, act{Q) = hold. 

As expected, well-formedness is preserved by structural equivalence: 
Proposition 2.1 (Well-formedness is up to =). Let P, Q systems s.t. P = Q, P is well-formed if it is Q. 

2.2 Type system 

In this subsection we present a type system for Bio/3 terms, which provides a formal procedure for check- 
ing well-formedness. As usual in type theory, terms are typed by judgements w.r.t. some environment. 

Definition 2.2 (Judgement). A type T C ACT is a finite set of mobility actions. 

A judgement for a Bio/3 term K (a system or a membrane) is of the form i"i ; T 2 h K : T, where the 
environment T\ ; T 2 is formed by two sets of names, such that Ti n T 2 = 0, each of which is written as a 
list of the form xo,...,x„. 

The two sets in the environment keep track of the free names in terms, and the number of times each 
name occurs: in our case, if n G Ti it occurs once, if n G T 2 then it occurs twice. Types keep track of 
mobility actions/co- actions occurring with free names. 

Moreover, let ok e and ok, be two predicates defined as follows: 

oMn^A^A,) ^ (r 1 ur 2 )n(A 1 UA 2 )=0 ok,(r,T,a) (z\ r ) ± = c\ r 

The first predicate checks if two environments share no names. The second predicate checks if the two 
types T, a pair-up correctly, i.e., they create correct actions-coactions pairs using names in T. 

The Bio/3 type system is shown in Figure 2. Intuitively, (empty) types empty systems and mem- 
branes, which have no names and no actions, (prot) deals with proteins (possibly) having free names 
which fit in the environment, and no actions. The rule (co-f) types the singleton membrane containing a 
co-fuse, which has a name x of rank 1 and the set containing itself as type, (action) is for the remaining 
actions: this rule checks that the subsystem contains no actions (thus enforcing action prefix of Defi- 
nition 2.1), adds afresh name x for the action, and create a type containing the action itself. The rules 
(v-prot) and (v-action) manage the restriction, in the first case v binds a name used by proteins, so the 
rule removes eventually x from names with rank 2 (remind that only names with or 2 occurrences can 
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be tied); in the second case there is a name x of rank 2 which connects a pair action-co-action, so the rule 
removes x from T2 as in previous case, and removes the actions from the type (such actions have no free 
names anymore). Last two rules are for composing systems and membranes, (par) puts side by side two 
subsystem, checking that the global name rank is no greater then 2 by ok e (enforcing graph-likeness) 
and that actions pairs with respective co-actions by ok t (action pairing); notice that common names (T) 
are promoted from rank 1 to 2. For (cell) is similar: we build a cell from subsystems and, beside the 
previous checkings, we impose that all names of rank 1 of P are contained in the ranked 1 names of S, 
i.e, T2 is empty; this enforces impermeability. 

Proposition 2.2 (Unicity of type). Let A" be a Bio/3 system or a membrane. If ri;]^ h K : x and 
Ai ; A2 h K : c, then T\ = Ai, T2 = A2 and x = c 

Theorem 2.1 (Well-formedness). A Bio/3 system P is well-formed if and only if Y\ ; T2 h P : X, for some 
environment Ti , T2 and type x. 

Proposition 2.3 (Subject congruence). Let P, Q be two Bio/3 systems and S, T two Bio/3 membranes, 
such that P = Q and S = T, then H ; T 2 h P : X if and only if Ti ; Y 2 H Q : X and Ai ; A 2 h S : o if and only 
if Ai;A2 h T : a, for some T2, Ai, A2 and x, G. 

From now on, we consider only well-formed terms, except when it is explicitly said. 

2.3 Semantics 

The operational semantics of a Bio/3 system is given by a reaction relation — ► on systems, defined by 
means of a rule system parametrized over basic protein reactions and mobility configurations. 

Defining a reaction semantics by means of a set of rules is the simplest way for a biologist to look at 
protein reactions, since it resembles usual chemical reactions. Moreover, the rule-based modeling leaves 
biologists a great degree of freedom in defining the behavior of a protein in solution. Such modeling 
freedom is somewhat limited by biochemical constraints. One of this was individuated and formalized 
by Danos and Laneve in [8]: the causality constraint. 

Biochemical reactions are either complexations (i.e. where low energy bonds are formed on two 
complementary sites) or decomplexations, possibly co-occurring with (de)activation of sites. Causality 
does not allow simultaneous complexations and decomplexations on the same site. This constraint is 
assured by the notion of ( anti- )monotonicity, which forces reactions not to decrease (increase) the level of 
connections of a solution [8]. Moreover, protein reactions should be closed on all well-formed membrane 
contexts, so that we are able to model reactions for proteins located in different membranes and cells. 

Let us now discuss about membrane level interactions, that is, rearrangement of the membrane struc- 
ture. As observed by biologists, these interactions result from complex nets of signaling pathways which 
induce a mechanical reshaping of biological membranes, allowing for example the formation of new 
vesicles, and hence new compartments. However, as noticed before, actual membrane transformations 
are limited and can be either a fuse or a pinch [6] ; what changes from situations to situations are the 
proteins involved in the signaling pathway, leading to the actual fuse or pinch. In order to formalize 
this situation we will distinguish between the effective membrane rearrangement, and what has trigger 
it. More precisely, a membrane reaction is split into two steps: a "preparation phase" where mobility 
actions (p or f) are introduced as prefixes in the system, and a "commitment phase" where the structural 
reconfiguration of nested compartments is actually performed. The commitment phase is formalized by 
rules defined once and for all by the framework, corresponding to the pinch and fuse processes of Fig- 
ure 1, but the preparation phase is specified by the modeler. In other words, one has only to describe 
when these interaction take place by indicating which are the suitable protein configurations that trigger 
a pinch or a fuse. In this way, each membrane interaction is given an explanation at the protein-level. 
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2.3.1 Protein reactions across multiple localities 

When dealing with protein reactions it is useful to think at complexes as single entities which can be 
affected in all their protein sub-units. An example of protein reaction we want to model is the following: 

D(\)*A(l+2 x )*lB(\ x + 2y)lC(\ y + 2)^^vz\D(\ z )*A(\ z + 2 x )*lB(\ x + 2y)lC(\ y + 2)^ 

A complexation reaction between protein D and the transmembrane ABC-complex has effect not only in 
the formation of a new protein bond between sites (D, 1) and (A, 1), but it toggles from hide to visible site 
(C, 2), which is not local to protein D (promotor of the reaction). This is a very common situation in pro- 
tein signaling transduction, indeed transmembrane receptors (as ABC-complex) propagate extracellular 
signals (as protein D) into the intracellular environment according to this mechanism. 
Let us now consider the following reaction: 

C(\ y + 2) * W + \A{l+2 x ) *D(1)H - Vz.(C(F + 2) * £5(1* + 2>') \A{V + 2 X ) *D(1 Z )^) 

This reaction differs from the previous one only for the direction in which the "signal" is propagated. 
There is no biological reason to distinguish between these two forms of reactions, indeed they are exactly 
the same protein reaction. In this sense, membranes have to be taken into consideration only in the way 
they partition solutions into distinct locations. 

In order to define protein interactions in nested systems we introduce (linear) contexts. Contexts 
generalize the Bio/3 syntax introducing variables, which are sorted into two categories, systems variables 
X and membrane variables Y, used as "holes" for systems and membrane, respectively. Sequences 
, . . . , K n are denoted by K, and concatenation by K, K' = K\ , . . . , K n , K[ , . . . , K' m . 

C[X;Y],D[X;Y] ::=o\A p (p) \X | p n ?5 | W[SlP^ | 

C[X>; Y'\ *D[X";Y"} | vn.C[X; Y] | lC[X; Y"]^ (System contexts) 

E [Y ] , F [Y ] : : = | A ap (p ) | p^-? 5 | | E [Y '] * F [Y>>] (Membrane contexts) 

where XeljG?, and let n, %' denote permutations, X = n(X',X"), Y = n'(Y', Y"). Context applica- 
tion C[P;S] is defined as expected, for P and 5 two sequences of systems and membranes, respectively. 

As in the case of the fc-calculus, protein reactions can be of two kinds: monotone and anti-monotone. 
Monotonicity is defined on a growing relation ([>) over Bio/3 protein solutions. Differently from K- 
calculus, the notion of "protein solution" in the Bio/3 Framework have to take into account system and 
membrane proteins, and the fact that they are in different locations. Thus, we define protein solutions as 
sequences of system groups and membrane groups of proteins, as follows: 

Definition 2.3 (Wide protein solutions). A group of system proteins is a Bio/3 system P of the form 
Ai(pi) * . . . *A n (p n ), where A, G 0P p (1 < i < n); similarly, a group of membrane proteins is a Bio/3 
membrane S of the form Bi((7i) *. . .-kB m (o m ), where Bj 6 &* a p (1 < J < m). 

A wide protein solution is a couple of sequences of system groups and membrane groups of proteins, 
which can be restricted on some names x, denoted as Yx.(P\, . . . | Si, . . . ,Sh). 

Informally, a group of proteins is a local protein solution, that is, all proteins in that group reside in 
the same location. Wide protein solutions are non-local and are formed as sequences of protein groups 
of two kinds according with the syntactic separation of membranes and systems. Actually, wide protein 
solutions are not Bio/3 terms, but we use them as specifications of how complexes are formed abstracting 
on localities. Often, we will shorten vx.(P \ S) as (P \ S) when \x\ = 0. Free names are defined as 
fn(vxi, . . .x n .(P | 5)) = fn(P,S) \{x u .. .x n }. 
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We are interested in connected wide protein solutions, that is, wide solutions where exists a path 
of bonds linking all the proteins in the solution. Closed connected wide protein solutions are actually 
protein complexes. Formally: 

Definition 2.4 (Connectedness). The connected wide protein solutions are defined inductively as: 

(P \s) <g | f) 

Vx.(P\S) fn((P\S))nfn((Q\f))^<d 



(A(p)|) <|fl(<r)> vxx.(P\S) (P,Q\S,f) 

(Pi,P,P 2 \S) (21) (P\S h S,S 2 ) (\T) 

fn((P h P,P 2 \S))nfn((Q\))^(d fn((P\S l ,S,S 2 ))nfn((\T))^(d 

(P l ,P*Q,P 2 \S) (P\S U S*T,S 2 ) 

A wide protein solution vx.(P \ S) is a wide complex if it is connected and fn(vx.(P \ S)) = 0. 
Definition 2.5 (Growing relation >). Let x be a set of fresh names. A growing relation > over interfaces 
is defined inductively as follows: 

iEi xn/«(p)=0 xhpoa xhp'oa' 
Jc h i > i xV- i\>i x h i\>i x x h p > p Jc h p + p' [> a + a' 
A growing relation \> over wide protein solutions is defined inductively as follows: 

xh (P\S)>{P' \S') xY- {Q\f)t>{Q \T) 

*h(o | )>(<>|) xh(|0)>(|0) x\- (P,Q\S,f )>(?>, Q' \S',T>) 

xh (P,P\S)>(Q,Q\f) ihp>0 x\- (P\S,S)\>(Q\f,T) xhp>a 

xh{P,P*A(p)\S)>(Q,Q*A(o)\f) xh{P\S,S*B(p))>(Q\f,T*B(o)) 

xh(P\S)>(Q,Q\f) fn{p)Qx xY~ (P \ S) D> (Q \ f ,T) fn(o)Qx 

xV- (P\S}> (Q,Q*A(p) | f) x\-(P\S)>(Q\ f,T*B(o)) 

Here, growing means possibly creating "new" bonds and/or proteins, and toggling sites from visible 
to hide or vice versa. Note that new bonds use names from x, that is, the following proposition holds: 
Proposition 2.4. If x h (P\S)>(Q\ f), then fn({Q | f))=fn({P \ S)) U x and x n fn({P \ S)) = 0. 

On the concept of growing we define monotonicity, which is very important to guarantee that the 
protein reactions respect the principle of causality. 

Definition 2.6 (Monotonicity and Anti-monotonicity). Let (P \ S) and (Q \ T) two wide protein so- 
lutions, we say that the pair ((P \ S),Vx.(Q \ T)) is monotone if x h (P \ S) > (Q \ T) and vx.(Q \ T) is 
connected. (vx.(Q \ T), (P \ S)) is anti-monotone if ((P \ S),vx.(Q \ T)) is monotone. 

We can now define what is a protein reaction specification, and the induced protein reactions. 
Definition 2.7 (Protein reaction specifications). A protein reaction specification is a set of monotone 
and/or anti-monotone pairs of wide protein solutions. 

Given a protein reaction specification the corresponding protein reactions are defined as follows: 

{(P\S),Vx.(Q\f))e<% (vx.{Q\f),{P\3))€& . t . 

— — — — '4-^ — (mon) v L -, — ^ — (anti-mon) 

C[P; S] vx.C[Q; T) vx.C[Q; T] C[P; S] 

where C[X;F] is a context such that vx.C[Q; T\ is well-formed (i.e., well-typed). 

In virtue of the generality of context C[X;F], we can adopt a simplified notation for reactions: (P \ 
S) — > vx.(Q | T) for monotone reactions and, similarly, vx.(Q \ T) —>■ (P \ S) for anti-monotone ones. 
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2.3.2 Membrane transport 

Membrane transformations are sorted into two groups: pinch (p) and fuse (f). They occur only with the 
interaction of an action t n and a co-action t„ (t G {p,f}). Pinch and fuse can be applied with different 
orientations, depending on where actions constructs are placed. Formally, the commitment rule are 

Vn%P*lvt%S*TlQ^lT\lSlPft*Qft (pinch-in) 

? 5* T I p n 5 P * Qtf - IS I Ptf * IT I Qtf (pinch-out) 

fn',lSlPtt*Un*nQti^lS*TlP* Qtt (fuse-hor) 

^ *Tlf n ;lSlPH*QH^P*lS*TlQH (fuse-ver) 

We can pinch inward or outward a sub-system depending on the direction where p and p- 1 are placed 
relatively to each other. Note that the continuations of p and p L are used to form the new compartment, 
either in (pinch-in) or (pinch-out). Fusions can be either horizontal (fuse-hor) or vertical (fuse-ver) 
according to the relative positions of compartments to be merged. Notice that, even in these reactions, f 
and f 1 - continuations are used to select which part of the system has to be merged. 

Mobility reactions respect bitonality in the sense of [4, 6], in fact the parity of nesting of P and Q is 
preserved in all these reactions, hence they preserve the bitonal coloring of those subsystems. 

In nature, membrane transport occurs when a certain protein configuration, allowing the mechanical 
reshaping of the membrane, is reached. Thus, in order to specify when a pinch or a fuse takes place, it is 
sufficient to specify which are the protein configurations leading to the "appearing" of mobility controls. 

Definition 2.8 (Membrane reaction specifications). A pinching configuration is a 5-tuple (P, P' ,S' ,S' ,Q) € 
IIxIlxIxExn, where IT, E are the sets of action-free systems and action-free membranes, respec- 
tively. A fusing configurations is 5-tuple (P,S,R,T,Q) 611x1x11x1x11. 

A membrane reaction specification is a pair (& p ,3£f) where £% p is a set of pinching configurations, 
and &f is a set of fusing configurations Given a membrane reaction specification (^ p ,^f), the corre- 
sponding mobility introduction reactions are defined as follows: 



(intro-p-in) 
(intro-p-out) 



(P,P',S,S',Q)e^ p fn(P)Cfn(S) fn(S) C fn(P,Q) n fresh 
P*P'*lS*S'*TlQ*R)^ Vn.(p n ° 9 P*P'*lp±° 9 S*S'*TlQ*R^) 

(P,P',S,S',Q)e% fn(P)Cfn(S) fn{S) Q fn{P,Q) n fresh 
Q*lS*S , *TlP*P , *Rft-> vn.(Q*lp^° 9 S*S'*Tlp n ° 9 P*P'*R^) 

(P,S,R,T,Q) £ gg f n fresh 
lS*S>lP*P>^*R*lT*T>lQ*Qi^^ vn.(f n °,lS*S>lP*pi^*R*lf ± *T*TilQ*Qi^) 

(intro-f-hor) 

(P,S,R,T,Q) €& f n fresh 
IT*T>IR*IS*S>IP*P>^*Q*Q>^^ vn.(lf^*T*T>lR*f n °,lS*S>lP*P>^*Q*Q>^) 

(intro-f-ver) 

Intuitively, the 5-tuples in & p and &f describe the systems configurations which trigger a mobility 
action. Notice that in the two pinch introduction rules, the first two pairs of terms (P,P') and (5,5') 
in the tuple correspond to two splits of sub-terms in the left hand side of the reaction rule: this is for 
distinguishing the part that is actually pinched (P and 5) from the part that just helps triggering reaction 
but do not change its position (P' and 5', which may be enzymes). The same holds for fuse introduction 
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rules, where P and S will be fused with the compartment where T and Q occur, whereas R just helps 
the fusion but will not change its position. In practice, actions are used for "freezing" the state of the 
subsystem involved in the reconfiguration, until the reconfiguration actually takes place. 

Informally, ^ p and ^ provide a template of suitable conditions, cause of which membrane recon- 
figuration can happen. Such conditions are to be defined by the modeler, who knows which are the 
requisites to make possible a pinch (either for vesiculation or endocytosis) or a fusion (either for fusion 
or endocytosis). This way of describing a model in the Bio/3 framework gives biologists (in general, 
modelers) the possibility to justify membrane reconfigurations at the protein-level, providing only the 
essential characteristics that a sub-system must satisfy to perform a mobility reaction. 

Summing up, to define a Bio/3 model, users have to provide the protein signature, the set of protein 
rules ignoring compartments, and two sets of configurations ^ p and M$ describing the situations 
for which a mobility reaction can occur. All other reactions are automatically provided by the frame- 
work. Hence, modelers can focus on protein-level aspects, leaving the framework to deal with biological 
compartments and membrane transport. 

Definition 2.9 (Bio/3 specification). A Bio/3 specification is a quadruple where & is a 

protein signature, M is a protein reaction specication, and (& p ,3!f) is a membrane reaction specication. 

Definition 2.10 (Bio/3 reactive system). Let ,&,& p ,&i) be a specification. The associated Bz'o/3 
reactive system is a pair (H act , — >), where H act is the set of well-formed Bio/3 systems over and — 
called the reaction relation, is the least binary relation over H act containing all protein reactions over 
mobility reactions pinch-in, pinch-out, fuse-hor, fuse-ver, and all introduction reactions over & p and <^f, 
and it is closed under structural equivalence, system composition, restriction of names, and compartment 
nesting. 

Closure of — > under system composition, restriction and compartment constructs, allow to focus on 
the actual "reacting parts" of the system. Clearly, — ► is not closed under action prefixes, because action 
prefixes "freeze" their continuation so reactions cannot happen in them. 

We can prove that for Bio/3 reactive systems subject reduction holds, that is, well-formedness for 
systems is preserved by reaction steps. 

Proposition 2.5 (Subject reduction). Let P, Q be two Bio/3 systems. If Ti;]^ h P : T and P — > Q, then 
ri;A2 h Q : a where either r2 = A2 and t = a, or r2 = A2,n and T = o + {t n ,t„} for some t € {p,f}. 

3 An example: Membrane Traffic 

In this section we give a formal Bio/3 description of how membrane traffic works. Here we take two 
important processes into account, which are fundamental in many transport mechanisms in cells: vesicle 
formation and vesicle docking. Vesicles form by budding from membranes. Each bud has distinctive 
coat protein on cytosol surface, which shapes the membrane to form a bubble. The bud captures the 
correct molecules for outward transport by a selective bond with a cargo receptor attached through the 
double layer to the coat protein complex, which is lost after the budding completes (see Figure 3(a)). 
Vesicle contents are released in the target membrane by fusing vesicle with it. 

To ensure that membrane traffic proceeds in an orderly way, transport vesicles must be highly selec- 
tive in recognizing the correct target membrane with which to fuse. Because of the diversity of membrane 
systems, a vesicle is likely to encounter many potential target membranes before it finds the correct one. 
Specificity in targeting is ensured displaying protein surface markers that identify them according to 
their origin and type of cargo (SNAREs). Vesicles are marked by v-SNAREs (vector SNAREs) which 
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Figure 3: (a) Vesicle formation, (b) Vesicle docking. (Pictures taken from [2]) 



selectively bind with some corresponding t-SNAREs (target SNAREs) located on target membrane (see 
Figure 3(b)). SNARE proteins have a central role both in providing specificity and in catalyzing the 
fusion of vesicles with the target membrane. 

Let & p = {R e i ,R c i ,C i ,AduCl,Sn e i ,Sn c i ,tSn^tSn'i} and 9> ap = {Rf ,Snf ,tSnf }, for i E {0, 1}. Trans- 
membrane proteins are represented by three sub-units in order to specify the direction they exhibit with 
respect to membranes into which are embedded. We denote with 7?- (l +2), ^"'(l + 2), 7??(1 +2) the 
three cargo receptor sub-units, respectively for cytosolic, trans-membrane, and extra-cytosolic; SNAREs 
(SnJ and tSn*) are defined in a similar way. Q(l) denotes a cargo molecule, Adj(l + 2) adaptin, and CI (I) 
clathrin. Proteins are sub-scripted with an index i £ {0, 1} to distinguish different families of proteins. 

We describe the BiojS model for membrane traffic giving the set of reaction rules for proteins and 
introduction of mobility. Protein reactions are defined as follows: 

(Q(l)*tf?(l+2*),/??(F + 2) |*S»(l* + 2>)) ^ vz.(Q(l z )*^(l z + 2- v ),^(l>' + 2) | /?f (l- v + 2>')) 
{R c i {\ x + 2)*Ad i {\ + 2) |) ^^vy.(R'i(\ x + 2?')*Adi(l y + 2) |) 
(Adi(\ x + 2)*Cl(l) | ) vy. (Adi(\ x + 2 y ) * Cl(V') |) 
vy.(R'i(l x + 2y)*Ad i (l y + r)*Cl(\ z ) |) (R'j(l x + 2)*Ad i (\+2 z )*Cl(l z ) |) 

(5^'(l' r + 2)*^(l+2 z ) | ) vy. (Sn c i ( l x + 2>> fcSn? {V + 2 Z ) \) 

Intuitively, (rec) describes the formation of a bond between the cargo molecule and the corresponding 
cargo receptor, which toggles to visible the hidden site of the receptor cytosolic sub-unit, (adpt) and 
(coat) allow for the creation of the clathrin coat which leads to the formation of the bud. (uncoat) frees 
the vesicle from the adaptin-clathrin coat. Finally, (snare) deals with the snare-mediated docking of the 
naked vesicle, by forming a bond between the v-SNARE and t-SNARE of the same family. Figure 3 
gives an informal idea of how protein reactions work. 

Introduction rules are defined by the two sets of pinching and fusing configurations: ^ p and Sfy . 

P = Q{l x ) *R e i (\ x + 2y)*Sn e i (V) P' = o 

<% P = {{P,P',S,S',Q)} where S = Rf{V + 2 w )*Sn? (l z + 2") S' = 

Q = R<!(r + T) *Adi(\ a + 2 b ) * Cl(l b ) * 5<(1" + 2) 

The configuration (P,P' ,S,S' \Q) allows the vesicle to detach when the adaptin-clathrin coat completely 
cover the forming bud. For sake of simplicity we consider coats of size one, i.e. with only one adaptin- 
clathrin complex (generalization to coat of size n is straightforward). 
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P = Q(\ x ) *R^\ x + 2 y )*Sn e i (\ a ) T = tSn? (l d + 2 e ) 

&f = {(P,S,R,T,Q)} where 5 = R? (\ y + 2 Z ) * Sn? (\ a + 2 b ) Q = tSn e i (\ e ) 

R = R c i (V + 2)* Sn'i{\ b + 2 C ) * tSn^Y + 2 d ) 

The configuration (P, S,R, T, Q) describe the situation when the vesicle is docked to the target membrane, 
hence the fusion can take place. 

One possible initial configuration can be defined as follows: 

Main = DOrg * Targo * Targ\ 

DOrg = Cl(l)*Cl(l)*Ad (l+2)*Ad l (l+2)*R c (l yo +2)*R'i(l yi +2)*Sn c Q (\ hQ + 2)*Sn\{\ hl +2)* 
IR™(\ X « + 2 ya )*R™(\ X0 +2 yi )*Sn™(\ aa + 2 ho )*Sn™(\^ +2 bl ) 

lRt ) (\+2 X0 )*R\(\+2 Xl )*Sn e (l a °)*Sn\(l a >)*C Q (l)*C l (\m 
Targ () = tSn c (l + 2 d °) * ^'(1* + 2 C ») ltSn e (l e °)^ 
Targx =tSn c l (\+2 do )*ltSn1\l dl +2 ei )ltSn e l (\ ei )^ 

In this case we have a single donor organelle which is able to pinch out two vesicles (using two different 
type of SNAREs). Such vesicle can approach the corresponding docking membrane. In Figure 3(b) is 
depicted a possible run of computation. 

4 Comparing with zc-calculus 

In this section we give a formal connection between Bio/3 and Danos and Laneve's fc-calculus [8]. This 
calculus has been taken as the formal model protein interactions in the definition of biobigraphs, the 
bigraphical model underlying the Bio/3 Framework. 
The syntax of fc-solutions is the following: 

S,T ::=0 | A(p) | S ,T | (x)(5) 

A ff-solution can be either the empty solution 0, or a protein A(p), or a group of solutions S , T, or a 
solution prefixed by name restriction (x) (S) with x G . 

As usual in process calculi, a structural equivalence (=) over processes is introduced, which is the 
least structural containing a-equivalence and satisfying the abelian monoid on composition and the scope 
extension and extrusion laws. Moreover we define the sets of free edge names fn(p) and fn(S), on 
interfaces and solutions, defined in the common way, where the new name constructor (x) (S) is the only 
name binder. We say a ff-solution S is closed iff fn(S) = 0. 

Connectedness of a ^-solutions is particular a simper case of the one seen before. 

Definition 4.1 (Connectedness and complexes). Connected K-solutions are defined inductively as: 

S ST fn(S)nfn(T)^<d S S = T 

Mp) Wsj f 

A ^-solution S is a complex if it is closed, graph-like, and connected. 

The semantic of fc-calculus is defined by means of a transition system built up by a user-defined set 
of monotone and anti-monotone rules & K . 
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Definition 4.2 (Growing relation). A growing relation [> over (("-solutions is defined inductively as: 

xhS\>T xhp>a xhS\>T fn(<j)Cx 
xhOoO xhS,A(p)i>r,A(a) xhS[>7\A(a) 

Again, monotonicity is a simplification of the one shown for Bio/3 . 

Definition 4.3 (Monotone reactions). Let L,R be two K"-solutions. L — ► (x)/? is a monotone reaction if 
x h L D> /?, where L, (Jc)/? graph-like and /? connected. (Jc)L — ► /? is antimonotone if /? — > (x)L is monotone. 

Finally, we are able to characterize a protein transition system: 

Definition 4.4 (Protein transition systems). Given a set 3% K of monotone and antimonotone reactions, 
we define a protein transition system (PTS) as a pair such that 5? is a set of JC-solutions, and ~> 

is the least binary relation over 5? such that 3% K C closed w.r.t. =, composition and name restriction. 

Now we are able to establish formally a connection between the Bio/3 framework and the ff-calculus, 
that is all protein reactions that can be take place in our systems are justifiably by means of protein 
reactions in a ff-solution. In order to do this, first we show a way to translate Bio/3 systems and reactions 
into k solutions and reactions ([—]). [[—J is defined as follows: 

[o] = o \A p { P )i=A p { P ) \p*Q\ = \p\AQ\ p™ = M,p>fl 

10] =0 \A ap {p)\=A ap {p) [S*T} = [Sl,[Tj [vn.P] = (n)([P]) 

lp n °,P} = IP} [f„ 5 P] = IP} U 5 S} = IS} [tf ] = 

Such encoding induces a type system over the fc-solution for checking if the solutions are graph-like. 
It can be given as follows: 

, v A£^ Vx£fn(p).\p,x\ <2 

lZ6rOJ 0;0hO {i£ fn{p) | |p,jc| = l};{jc€/n(p) | |p,jc| = 2} h A(p) lPr ° j 

, , r i; r 2 h5 x^n r!,r ; r 2 h5 A 1 ,r ; A 2 hr (r 1 ur 2 )n(A 1 uA 2 ) =0 

(res) — — (par) 

1 ; r i; r 2 \{x}h(x)5 FuA^M^h S,T {W > 

The following properties stating a statical relation between Bio/3 and K hold. 

Proposition 4.1 (Syntax). 1. For S a fc-solution, if Ti;r 2 hS and Ai;A 2 h 5thenTi=Ai andT 2 =A 2 . 

2. For 5 a fc-solution, S is graph-like iff Ti ; T 2 h 5 for some Ti , T 2 . 

3. For S, T two K-solutions, if 5 = T, then Ti ; T 2 h 5 iff Ti ; T 2 h T . 

4. For P a Bio/3 system, if Ti ; T 2 h P : t, then Ti ; T 2 h [P] . 

Finally, we can state and prove the following semantics/dynamic relation between Bio/3 and k: 
Proposition 4.2 (Semantics). 1 . For S, T two fc-solutions, if Ti ; T 2 h 5 and S ~> 7\ then H ; T 2 h P. 

2. For a wide monotone protein reaction ((P,S), Vx.(P',S')) 6 and for every C[X,?] such that 
C[P,S] is well-formed, then [C[P,5]] ~» [vx.C[P', 5']]. 

3. For a wide anti-monotone protein reaction (vx.(P,S), (P',5')) G for every C[X,?] such that 
Vx.C[P,S] is well-formed, then [vx.C[P,5|] ~> [C[P',5']]. 

Intuitively, this proposition states that every protein transition performed in the Bio/3 framework 
is justified by a corresponding transition in the fc-calculus. The converse does not hold, due to the 
compartment confinements: not all reactions of the translation of a Bio/3 system which are possible in K- 
calculus, correspond effectively to some wide (anti)monotone protein reactions in the Bio/3 framework. 
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5 Related works 

In literature several other calculi, covering both protein reactions and membrane interactions, have been 
proposed. Two particularly interesting cases are the biofc-calculus [9] and the "^-calculus [7]. In this sec- 
tion we compare our framework with these calculi analyzing analogies and differences of the approaches. 

biofc-calculus can be seen as a version of the ff-calculus extended with compartments. Compared 
to Bio/3, the most remarkable similarities are at the syntactic level. Both calculi have a compartment 
constructor which embeds proteins into membranes, and similar well-formedness conditions to restrict 
to only biological relevant processes. Both calculi are sound at protein level in the sense of the ff-calculus. 
In both cases preservation of well-formedness during reactions is proved, although biofc-calculus does 
not have a formal type system as Bio/3. 

Several differences emerge at the semantic level, though. An important design choice distingushing 
bioff and Bio/3 concerns how the behaviour of a compartment is specified. In bioff, the mobility capa- 
bilites of a compartment are indicated by its name, which appears explicitly in the rearrangement rules 
which can be applied. The result of a membrane interaction is one or more new membranes whose names 
have to be specified in the rearrangement rule. This means that the rule defines also the future (mobility) 
behaviour of the resulting systems; hence it is up to the modeller to specify correctly the evolutions of 
names of membranes. 

On the other hand, in Bio/3 compartments have no names: membranes are just containers and their 
behaviour is completely determined by the interactions of proteins floating in the aqueous solutions and 
lipid bi-layers. Therefore, a modeller has only to provide protein-level rules and define the configurations 
which trigger a mobility reaction. One does not have to declare how a membrane will evolve (i.e., "which 
is the next name" in biofc), because the behaviour is defined by the rules of the framework. 

Due to these radically different design choices, a formal comparison at the semantic level between 
biofc-calculus and Bio/3 is not easy. Moreover, the bioff semantics is given by means of an LTS, instead 
Bio/3 has a reaction (reduction) relation. A possible future work is to derive an LTS with associated 
congruential bisimilarity, taking advantage of the theory of bigraphs. 

"^-calculus Both 'rf -calculus and Bio/3 frameworks are founded on bigraphical reactive systems, which 
does not only confer a formal graphical representation but provides also a powerful categorical theory 
for automatically constructing labelled transition systems (LTSs) whose bisimilarities are always con- 
gruences. Even though ^-calculus and Bio/3 share common sources of inspiration, they differs in many 
aspects, one of which is the way compartment reconfiguration are performed. 

^-calculus adopts meta-biological "gates" to represent a semi-fused intermediate state in membrane 
fusion or fission. Gates act like channels, allowing protein diffusion between compartments, and it does 
not allow for phagocytosis. Moreover, one of the aim of Bio/3 framework is to characterize higher-level 
membrane interactions with lower-level protein reactions, and to be adequate with respect to the K- 
calculus; a formal connection between ^-calculus and ff-calculus has not been investigated yet. Finally, 
in the ^-calculus, mobility and protein reactions are more intertwined, since membrane interactions are 
in many steps: introduction of gates, some protein diffusions, and deletion of gates. This mechanism does 
not provide an explicit, logical separation between compartment reconfiguration and protein reactions. 
On the other hand, the Bio/3 framework aims to an explicit distinction between protein reactions and 
mobility actions, as implemented by means of introduction rules, which freeze the subsystems involved 
in the rearrangement until the structural reconfiguration is actually performed. 
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6 Conclusion 

In this paper we have presented the Bio/3 Framework, a language for both (lower-level) protein and 
(higher-level) membrane interactions of living cells. A typed language has been given for describing 
only biologically meaningful compartments and proteins. 

In the Bio/3 framework membrane reconfigurations are logically different from protein reactions, 
but the connection between the two aspects is formally specified. In fact, membrane activities can be 
motivated by a formal biological justication in terms of protein interactions. Notably, protein interactions 
are adequate with respect to the ff-calculus. 

As an example application of this framework, we have modelled the mechanism of membrane traffic, 
which consists in formation, transportation and docking of vesicles. 

Related and future work. Due to lack of space, we have omitted a formal connection between the pro- 
posed framework and its inspiring bigraphical model. Another interesting aspect to investigate is the 
bisimilarity definable via the IPO construction on bigraphs [10], which is always a congruence. This can 
be useful e.g. in synthetic biology, for instance, to verify if a synthetic protein (or even a subsystem to be 
implanted) behaves as the natural one. We plan to describe such results in a forthcoming work. 

Another interesting investigation is to give a formal comparison of Bio/3 with the other two frame- 
work present in literature, such as biofc and "if -calculus. In particular, it could be useful to understand the 
differences of ^-calculus rule generators and our introduction rules, even because they are both inspired 
by bigraphical wide reactions. 

As a next step, a stochastic version of the Bio/3 framework should be analyzed for simulation pur- 
poses. One way could be trying to adapt the extended Gillespie simulation algorithm [11], that works for 
calculi with multi-compartments. 
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